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Abstract. Borg-type uniqueness theorems for matrix-valued Jacobi opera- 
tors H and supersymmetric Dirac difference operators D are proved. More 
precisely, assuming reflectionless matrix coefficients A, B in the self-adjoint 
Jacobi operator H = AS + + A~ S~ + B (with the right/left shift opera- 
tors on the lattice Z) and the spectrum of H to be a compact interval [E_ , E+] , 
E— < -E+, we prove that A and B are certain multiples of the identity matrix. 
An analogous result which, however, displays a certain novel nonuniqueness 
feature, is proved for supersymmetric self-adjoint Dirac difference operators D 
with spectrum given by [ - E 1 ^ 2 , -E\[ 2 \ U [E 1 ^ 2 , E 1 ^ 2 ] , < £_ < E+. 

Our approach is based on trace formulas and matrix-valued (exponential) 
Herglotz representation theorems. As a by-product of our techniques we obtain 
the extension of Flaschka's Borg-type result for periodic scalar Jacobi operators 
to the class of reflectionless matrix- valued Jacobi operators. 



1. Introduction 

As discussed in detail in [23], while various aspects of inverse spectral theory 
for scalar Schrodinger, Jacobi, and Dirac-type operators, and more generally, for 
2x2 Hamiltonian systems, are well-understood by now, the corresponding theory 
for such operators and Hamiltonian systems with m x m, m £ N, matrix-valued 
coefficients is still largely a wide open field. A particular inverse spectral theory 
aspect we have in mind is that of determining isospectral sets (manifolds) of such 
systems. In this context it may, perhaps, come as a surprise that even determining 
the isospectral set of Hamiltonian systems with matrix-valued periodic coefficients 
is an open problem. The present paper is a modest attempt toward a closer in- 
vestigation of inverse spectral problems in connection with Borg-type uniqueness 
theorems for matrix-valued Jacobi operators and Dirac-type finite difference sys- 
tems. It should be mentioned that these problems are not just of interest in a 
spectral theoretic context, but due to their implications for other areas such as 
completely integrable systems (e.g., the nonabelian Toda and Kac-van Moerbeke 
hierarchies), are of interest to a larger audience. 

Before we describe the content of this paper in more detail, we briefly comment on 
background literature for matrix-valued Jacobi and Dirac-type difference operators. 
Spectral and Weyl-Titchmarsh theory for Jacobi operators can be found in [4, Sect. 
VII.2], [20], [47], [49, Ch. 10] and the literature therein. The case of Dirac finite 
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difference operators was discussed in detail in [12]. Deficiency indices of matrix- 
valued Jacobi operators are studied in [35]-[37]. Inverse spectral and scattering 
theory for matrix-valued finite difference systems and its intimate connection to 
matrix- valued orthogonal polynomials and the moment problem are treated in [1], 
[2], [4, Sect. VII.2], [16]-[18], [21], [39], [40], [45], [46], [47], [49, Ch. 8], [50]. A 
number of uniqueness theorems for matrix-valued Jacobi operators were proved 
in [23]. Finally, connections with nonabelian completely integrable systems are 
discussed in [5], [6], [44], [48], [49, Chs. 9, 10]. 

Let C(Z) rxs be the space of sequences of complex r x s matrices, r, s G N. We 
also denote No = N U {0} in the following. 

The matrix- valued Jacobi operators H in £ 2 (Z) m discussed in this paper are of 
the form 

H = AS + +A~S~ +B, V(H)=f(Z) m . (1.1) 
Here S 1 * denote the shift operators acting upon C(Z) rxs , r,seN, as 

S ± f(-) = f ± (-) = f(- ±1), feC(Z) r *', (1.2) 

and A = {A(k)} kez G C(Z) mxm , A(k) > for all k G Z, B = {B(k)} keZ G 
C(Z) mxm , B(k) = B(k)* for all k G Z and m G N. 

The Dirac-type finite difference operators D in £ 2 (Z) m © £ 2 (Z) ro studied in this 
paper are of the form 

D = S„ + X=( _ ° pS+ n X *), V(D)=e 2 (Z) m (B£ 2 (Z) m , (1.3) 



K p-s- + x o 

where p = {p(k)} keZ C C(Z) mxm and S p and X arc of the block form 

P S+\ (0 X 



s ^{p-s- J' x ={x o) ^ 

with p(k),x(k) G C mxm invertible for all k G Z. Moreover, following [12], we may 
assume without loss of generality that for all k G Z, p(fc) is a positive definite 
diagonal mxm matrix (cf. Remark 5.3 for details)). 

Additional assumptions on the coefficients A(k), B(k), p(k), and x(fc), k 6 Z, 
will be formulated in Sections 2 and 5, respectively. 

The difference operators H and D represent the natural matrix-valued general- 
izations of Lax operators arising in connection with Kac-van Moerbeke and Toda 
lattices (cf. [8] , [22] and the references therein) and hence lead to nonabelian Toda 
and Kac-van Moerbeke hierarchies of completely integrable nonlinear evolution 
equations. 

Next we briefly describe the history of Borg-type theorems relevant to this paper. 
In 1946, Borg [7] proved, among a variety of other inverse spectral theorems, the 
following result. (We denote by a(-) and <7 e ss( - ) the spectrum and essential spectrum 
of a densely defined closed linear operator.) 

Theorem 1.1 ([7]). Let q G L 1 1 oc (R) be real-valued and periodic. Leth = —-£rz+q be 
the associated self-adjoint Schrddinger operator in L 2 (R) and suppose that a(h) = 
[eo, oo ) for some eo G M. Then q is of the form, 

q(x) = e for a.e. i£l. (1.5) 

Remark 1.2. Traditionally, uniqueness results such as Theorem 1.1 are called 
Borg-type theorems. This terminology, although generally accepted, is a bit un- 
fortunate as the same term is also used for other theorems Borg proved in his 
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celebrated 1946 paper [7]. Indeed, inverse spectral results on finite intervals in 
which the potential coefficient (s) are recovered from several spectra were also pio- 
neered by Borg in [7] and theorems of this kind are now also described as Borg-type 
theorems in the literature, see, e.g., [41]-[43]. 

A closer examination of the proof of Theorem 1.1 in [13] shows that periodicity 
of q is not the point for the uniqueness result (1.5). The key ingredient (besides 
a(h) = [eo, oo) and q real-valued) is the fact that 

for all x G M, £(A,x) = 1/2 for a.c. A G a ess (h). (1.6) 

Here £(A, x), the argument of the boundary value ^(A+zO, x) of the diagonal Green's 
function of h on the real axis (where g(z,x) — (h — zl)^ 1 (x,x)), z G C\cr(h)), is 
defined by 

£(A, x) = 7T _1 limIm(ln(g(A + ie, x)) for a.e. A G R. (1.7) 

Real- valued periodic potentials are known to satisfy (1.6), but so do certain 
classes of real- valued quasi-periodic and almost-periodic potentials q. In particular, 
the class of real- valued algebro-geometric finite-gap KdV potentials q (a subclass of 
the set of real- valued quasi-periodic potentials) is a prime example satisfying (1.6) 
without necessarily being periodic. Traditionally, potentials q satisfying (1.6) are 
called reflectionless (see [11], [13] and the references therein). 

The extension of Borg's Theorem 1.1 to periodic matrix- valued Schrodingcr oper- 
ators was proved by Depres [15]. A new strategy of the proof based on exponential 
Hcrglotz representations and a trace formula for such potentials, as well as the 
extension to reflectionless matrix- valued potentials, was obtained in [13]. 

The direct analog of Borg's Theorem 1.1 for scalar (i.e., m = 1) periodic Jacobi 
operators was proved by Flaschka [19] in 1975. 

Theorem 1.3 ([19]). Suppose a and b are periodic real-valued sequences in £°°(Z) 
with the same period and a{k) > 0, k G Z. Let h = aS + +a~ S~ +b be the associated 
self-adjoint Jacobi operator in £ 2 (Z) and suppose that a(h) = [E_,E + ] for some 
E- < E + . Then a and b are of the form, 

a(k) = {E+-E_)/4, b{k) = (E-+E+)/2, k G Z. (1.8) 

The extension of Theorem 1.3 to reflectionless scalar Jacobi operators is due to 
Teschl [53, Corollary 6.3] (see also [54, Corollary 8.6]). As one of the principal 
results in this paper we will extend Theorem 1.3 to matrix- valued reflectionless 
Jacobi operators H of the type (1.1) in Section 4. This extension (and the special 
case of periodic matrix-valued Jacobi operators) is new. In addition, we will prove 
a Borg-type theorem for the supersymmetric Dirac difference operator D in (1.3) 
which is new even in the simplest case m = 1. The latter displays an interesting 
nonuniqueness feature which has not previously been encountered with Borg-type 
theorems. 

In Section 2 we review the basic Weyl-Titchmarsh theory and the corresponding 
Green's matrices for matrix-valued Jacobi operators on Z and on a half-lattice. 
Section 3 is devoted to asymptotic expansions of Weyl-Titchmarsh and Green's 
matrices as the (complex) spectral parameter tends to infinity. Section 4 contains 
the derivation of a trace formula for Jacobi operators and one of our principal 
new results, the proof of a Borg-type theorem for matrix-valued Jacobi operators. 
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Finally, Section 5 presents a Borg-type theorem for supersymmetric Dirac-type 
difference operators. 

2. Weyl-Titchmarsh and Green's Matrices for Matrix- Valued 

Jacobi Operators 

In this section we consider Weyl-Titchmarsh and Green's matrices for matrix- 
valued Jacobi operators on Z and on a half-lattice. 

We closely follow the treatment of matrix- valued Jacobi operators in [23]. As 
the basic hypothesis in this section we adopt the following set of asumptions. 

Hypothesis 2.1. Let in G N and consider the sequences of self-adjoint mxm 
matrices 

A = {A(k)} keZ g C(Z) mxm , A(k) > 0, k e Z, 

B = {B(k)} keZ eC(Z) mxm , B(k) = B(k)\ keZ. (2.1) 

Moreover, assume that A(k) and B(k) are uniformly bounded with respect to k £ Z, 
that is, there exists a C > 0, such that 

||A(fc)|| C mxm + \\B(k)\\ C mxm < C, k £ Z. (2.2) 

Next, denote by S ± the shift operators in C(Z) rxs , r, s G N, 

(S ± g)(k) = g ± (k) = g(k ± 1), g G C(Z) r * s , k G Z. (2.3) 

Given Hypothesis 2.1, the matrix-valued self-adjoint Jacobi operator H with do- 
main £ 2 (Z) m is then defined by 

H = AS + +A~S~ +B, V(H)=l 2 (Z) m . (2.4) 

Because of hypothesis (2.2), H is a bounded symmetric operator and hence self- 
adjoint. In particular, the difference expression AS + + A~ S~ + B induced by (2.4) 
is in the limit point case at ±oo. We chose to adopt (2.2) for simplicity only. Our 
formalism extends to unbounded Jacobi operators and to cases where the difference 
expression associated with (2.4) is in the limit circle case at +oo and/or — oo (cf. 
[12]). We just note in passing that without assuming (2.2), the difference expression 
AS + + A~S~ + B is in the limit point case at ±oo if X^fc^fe ll^-(^)llc"» ><m = 00 ( see ' 
e.g., [4, Theorem VII. 2. 9]). Whenever we need to stress the dependence of H on 
A, B we will write H(A, B) instead of H . 

The Green's matrix associated with H will be denoted by G(z, k,£) in the fol- 
lowing, 

G(z,k,e) = (H-zI m )- 1 (k,£), zeC\R,k,£eZ. (2.5) 

Next, fix a site k E Z and define mxm matrix-valued solutions <f>(z, k, ko) and 
9(z, k, ko) of the equation 

A(k)tp(z, k+l)+A(k-l)tp{z, k-l)+{B(k)-zI m )tp{z, k) — 0, z G C, k e Z, (2.6) 

satisfying the initial conditions 

6(z,k ,k ) = <f)(z,k + l,fc ) = I m , <f)(z,k ,k ) = 0(z,k o + l,k ) = 0. (2.7) 

One then introduces mxm matrix- valued Weyl-Titchmarsh solutions ip±(z, k, ko) 
associated with H defined by 

ip±(z,k,k ) = 6(z,k,ko) - cj>{z,k,k )A{k y 1 M±(z,ko), z G C\R, k G Z, (2.8) 
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with the properties 

A(k)ip±{z, k + 1, fa) + A(k - l)V±(-2, k-1, ko) + (B(k) - zl m )^±(z, k, k ) = 0, 

z E C\R, k G Z, (2.9) 

^±(z r ,k Q ) e£ 2 {{k ,±oo)r)Z) mxm , zeC\R, (2.10) 



where M±(z, ko) denote the half- line Weyl-Titchmarsh matrices associated with H. 
Since by assumption AS + + A~ S~ + B is in the limit point case at ±00, M±(z, ko) 
in (2.8) arc uniquely determined by the requirement (2.10). We also note that by 
a standard argument, 

det((f)(z, k, k )) + for all k G Z\{fc } and z E C\R, (2.11) 

since otherwise one could construct a Dirichlet-type eigenvalue z € C\R for H 
restricted to the finite segment ko + 1, . . . , k — 1 for k > ko + 1 of Z (and similarly 
for k < ko — 1). Thus, introducing 

M N {z,k ) = -^(z.JV.fcoJ-^^.JV.fco), z€C\l, iVeZ\{M, (2.12) 

one can then compute M±(z, ko) by the limiting relation 

M±(z,k )= lim M N (z,ko), z G C\M, (2.13) 

TV— »±oo 

the limit being unique since + A~S~ + B is in the limit point case at ±00. 
Alternatively, (2.8) yields 

M±{z,k ) = -A(k )iP±{z,k + l,k ), zeC\I. (2.14) 

More generally, recalling 

det{i>±(z,k,ko)) + for all k G Z and z G C\M, (2.15) 

by an argument analogous to that following (2.11), we now introduce 

M±(z,k) = -A(k)i>±{z,k + l,k )il>±{z,k,ko)- 1 , z G C\R, fc G Z. (2.16) 

One easily verifies that M±(z,k) represent the Weyl-Titchmarsh M-matrices as- 
sociated with the reference point k G Z. Moreover, one obtains the Riccati-type 
equation for M±(z,k), 

M ± {z,k) + A(k-l)M ± (z 1 k-l)- 1 A(k-l) + zI m -B{k)=0, z e C\R, fcGZ 

(2.17) 

as a result of (2.9). For later reference we summarize the principal results on 
M±(z,k ) in the following theorem. (We denote as usual Rc(M) = (M + M*)/2, 
Im(M) = (M* - M)/(2i), etc., for square matrices M.) 

Theorem 2.2 ([3], [9], [30]- [33], [38]). Assume Hypothesis 2.1 and suppose that 
z G C\WL, and ko G Z. Then, 

(i) ±M±(z, ko) is a matrix-valued Herglotz function of maximal rank. In particular, 

Im(±M±(z, ko)) > 0, ze C+, (2.18) 

M±(f,feo) = M±(z,ko)*, (2.19) 

rank(M ± (z, fe )) = m, (2.20) 

limM ± (A + ie, fc ) exists for a.e. A G M. (2.21) 
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Isolated poles of ±M±(z,ko) and ^fM±(z,ko)~ 1 are at most of first order, are real, 

and have a nonpositive residue. 

(ii) ±M±(z,k ) admit the representations 

±M±(z,ko) = K±(k ) + L ± (k )z + £dQ±(\,k ) ^-1_- T A^ (2.22) 



C±(k )+ / d\Z±(\,k ) 



A 



A - z 1 + A 2 



(2.23) 



= exp 
where 

K±(ko) = K ± (k y, L + (k ) = o, Mfc ) = Mfc )*, / ll ^ 1 ± ? A ; 2 fco)l1 < oo, 

(2.24) 

C±(ko) =C ± (koT, O^S ± (-,fc )^/ m a.e. (2.25) 
Moreover, 

Q±((X,n],ko) = limlim- / dvIm(±M±(v + is, k )), (2.26) 

610 ej.0 7T J x+S 

S±(A,fc ) =lim-Im(ln(±M ± (A + i£,fc )) for a.e. Ael. (2.27) 

ej.0 7T 

Next, we define the self-adjoint half-line Jacobi operators H±^ ko on £ 2 ([ko, ±oo) n 
Z) m by 

#±,fc = P±,k (> HP± M \ e2({kQ ±oc)nZ)m , k G Z, (2.28) 

where P±.k are the orthogonal projections onto the subspaces £ 2 ([k , ±oo) n Z) m . 
In addition, Dirichlet boundary conditions at k =p 1 are associated with -H±,fc ; 

(H +M f)(k ) = A(fco)/+(fc ) + B(k )f(k ), 

(H- M f)(k ) = A-(k )f-(k ) + B(k )f(k ), 

(H±, k J)(k) = A(k)f + (k)+A-(k)f-(k)+B(k)f(k), fc|fc ±l, (2.29) 
/e»(%„) = / 2 (N,±™)nz) ra 

(i.e., formally, /(fc T 1) = 0). We also introduce ra-functions m±(z, k n ) associated 
with i?±,fc by 

m±(z, fc ) = Qk (H±,k - zI m y 1 Qk , (2.30) 
= G±, fco (z,A:o,fco), z£C\R, fc eZ. (2.31) 

Here Qfc are orthogonal projections onto the m-dimensional subspaces £ 2 ({ko}) m , 
fc 6 Z and 

G±, fco (z,M) = (#±,fc -^^(M), zeC\M, fc,£eZn[fc ,±oo) (2.32) 

represent the Green's matrices of H±,k„. 

In order to find the connection between m±(z,ko) and M±(z,ko) we briefly 
discuss the Green's matrices G±^ {z, k, £) and G(z, k, £) associated with £f±,fc and 
H next. 

First we recall the definition of the Wronskian W(f, g){k) of two sequences of 
matrices f(-),g(-) G C(Z) mxm given by 

W(f,g)(k) = f(k)A(k)g(k + l)-f(k + l)A(k)g(k), k G Z. (2.33) 
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We note that for any two matrix-valued solutions <p(z, •) and ip(z, •) of (2.6) the 
Wronskian W(<p(z, •)*, ip(z, ■))(&) is independent of k G Z. 

In complete analogy to the scalar Jacobi case (i.e., m = 1) one verifies, 

q i-i> + (z,k,k -l)A(k -l)-^(z,£,k -ir, £<k, 

+ ' {-^k^o-^Aiko-iy^+izJ^o-l)*, £>k, { ' ' 

z G C\R, k,l G Zn [fc ,oo). 

G_ fz fc fl = <! ^( z ' fc ' /c o + l)A(fc + l)-V-(^^,fco + l)*, , . 

_ ^ olZ ' ' ' \V-(«,*,*0 + l)^(*o + l)~V(^^*0 + l)*. ^>fc, 

z G C\R, G Zn(-oo,fc ]. 

Similarly, using the fact that 

^ + (z,k,k )[M_(z,k ) - M+^^o)]" 1 ^-^^,^)* 
= V-(-2,A;,fco)[M_(z,fco) - M+(z,fc )rV+(z,£;,fc )* 
= [M_(z, fc) - M + (z, fe)] _1 , fcGZ, (2.36) 

and 

fc + 1, fco)[M_(z, fc ) - M + (z, fco)]-V-(«, k, k )* (2.37) 
- i4(fe)^_(«, k + 1, fc )[M_(z, fc ) - M+(z, fco)] _1 V+(^ fc, M* = / m , fee Z, 
one verifies 

G , k t) = U+(z,k,k )[M_(z,k )-M + (z,k )}- 1 i>-(z,e,k y, e<k, 

\^(z,k,k )[M^(z,k )-M + (z,k )}- 1 ij + (z,£,k r, £>k, 

z G C\R, k,£eZ. 

Using (2.7), (2.9), (2.17), and (2.31), one infers that the Weyl-Titchmarsh ma- 
trices M±(z,k) introduced by (2.13) (resp. (2.14)) and the m-functions m±(z,k) 
defined in (2.30) are related by 

M+(z,k) = -m+(z,/c)" 1 - zl m + B(k), z G C\M, k G Z (2.39) 

and 

M-(z,k) = TO_(z,fc)-\ zgC\R, keZ. (2.40) 
In analogy to (2.17), m±(z, k) also satisfy Riccati-type equations of the form 
A(k - l)m+(z,k)A(k - l)m+{z,k - 1) + (zl m - B(k - l))m+{z, k - 1) + I m = 0, 

z G C\R, k G Z (2.41) 

and 

A(k - l)m-(z, k - l)A(k - l)m_(z, k) + (zl m - B(k))m,-(z, k) + I m = 0, 

z G C\R, k G Z. (2.42) 

Next, we introduce the 2m x 2m Weyl-Titchmarsh matrix M(z, fc) associated 
with the Jacobi operator H in £ 2 (Z) m by 

M(z, fc) = (M jtj , (z, fc)) ^ jl=1 2 , z G C\R, keZ, (2.43) 
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where 

M hl (z,k) = [M-(z,k) -M + (z,fc)]~ 1 , (2.44) 

Mi, 2 («, k) = 2- x [M_(z, fc) - M+(z, k)\- x [M-{z, k) + M+(z, k)], (2.45) 

M 2 ,i(z, k) = 2- x [M_(z, k) + M+(z, k)][M_(z, k) - M+(z, fc)]~\ (2.46) 

M 2 , 2 (z, fc) = M±(z, fc)[M_(«, k) - M+(z, fc)]- 1 M T (z, fc). (2.47) 

One verifies that M(z, k) is a 2m x 2m Herglotz matrix with the following prop- 
erties: 

Theorem 2.3 ([3], [9], [30]-[33]). Assume Hypothesis 2.1, z G C\R, and k G Z. 
T/ien, M(z, k ) is a matrix-valued Herglotz function of rank 2m with representations 

M(z, ko) = K(k ) + J dfi(A, feo) - ) ( 2 ' 48 ) 



exp 



C(fc )+ f dXT(X,k )(-^- 



A 



A - z 1 + A 2 



(2.49) 



where 



Moreover, 



dSl(X,k ) 

TTa 5 

C(fco) - C(fco)*, < T( • , k ) ^ I 2m a.e. (2.51) 



K(k ) = K(k y, I <oo, (2.50) 



fi((A,/x],A;o) = limlim- / di/ Im(M(i/ + ie, fc )), (2.52) 

<510 ej.0 7T J A+(5 

T(A,fc ) =lim-Im(ln(M(A + i£,fc ))) for a.e. X G R. (2.53) 

ej.0 7T 

Remark 2.4. We note that the Weyl-Titchmarsh matrix M(z, fc) is related to the 
Green's matrix associated with the Jacobi operator H by 

ft j (fc) W.*>ft J^HGl 7 ^)' (2 - 54) 



where 



With M(z,k) defined in (2.43), the following uniqueness theorem holds. 

Theorem 2.5. Assume Hypothesis 2.1 and let ko G Z. T/ien £/ie 2m x 2m Weyl- 
Titchmarsh matrix M(z, k ) for all z G C + uniquely determines the Jacobi operator 
H and hence A = {A(k)}k^z and B = {B(k)}kez- 

Perhaps the simplest way to prove Theorem 2.5 is to reduce it to knowledge of 
M±(z, k ), and hence by (2.39) and (2.40) to that of m±(z, k ) for all z G C+. (Wc 
note that the knowledge of B(k ), which is required according to (2.39), can be 
determined from the asymptotics of M_(z, ko) in (3.8).) Using the standard con- 
struction of orthogonal matrix-valued polynomials with respect to the normalized 
measure in the Herglotz representation of m±(z, ko), 



m±(z,k ) = / dv±(X,k ) (A - z) \ z G C+, / dv±{\, k ) = I m , (2 
Jr Jr 



56) 
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allows one to reconstruct A(k), B(k), k 6 [ho, ±oo)nZ from the measures dv±(X, ko) 
(cf., e.g., [4, Section VII.2.8]). More precisely, 

A(k Q ±k) = / AP ±;fc (A,fco)^±(A,fco)P±,fe + i(A,fc )*, 

B(k ±k) = [ AP ±i fe(A,fco)^±(A,fco)P±,fe(A,/e )*, keN , (2.57) 

where {-P±,fc(A, fco)}fceN is arL orthonormal system of matrix-valued polynomials 
with respect to the spectral measure dv±(\, k ), with P±,k(z, k ) of degree k in z, 
P±.o(z,k ) = I m - One verifies, 

P +ife (z,fco) = <t>(z,k + k,k -l), (2.58) 
P_ :fe (z,fco) = 6{z,k - k,k ), k g {-1}UN , fc G Z, z G C, (2.59) 

with (j)(z,k,k ) and 9(z,k,k ) defined in (2.6), (2.7). 

Given these preliminaries and introducing the diagonal Green's matrix by 

g(z,k) = G(z,k,k), z G C\R, k e Z, (2.60) 

we can also formulate the following uniqueness result for Jacobi operators obtained 
in [23]. 

Theorem 2.6 ([23]). Assume Hypothesis 2.1 and let k E Z. Then any of the 
following three sets of data 

(i) g{z, ko) and G{z, ko, ko + 1) for all z G C+; 

(ii) g{z, ko) and [G(z, ko, ko + 1) + G(z, ko + 1, ko)] for all z G C+; 
(m) (/(z, fco), g(z, ko + 1) for all z G C+ and j4(fco); 

uniquely determines the matrix-valued Jacobi operator H and hence A = {A(k)}kei, 
and B = {B(k)} keZ . 

The special scalar case m = 1 of Theorem 2.6 is known and has been derived in 
[53] (see also [54, Sect. 2.7]. Condition (iii) in Theorem 2.6 is specific to the Jacobi 
case. In the corresponding Schrodinger case, the corresponding set of data does 
not even uniquely determine the potential in the scalar case m = 1 (see, e.g., [27], 
[53]). 



3. Asymptotic Expansions of Weyl-Titchmarsh and Green's Matrices 

In this section we use Riccati-type equations to systematically determine norm 
convergent expansions of Weyl-Titchmarsh and Green's matrices as the spectral 
parameter tends to infinity. 

We start again with the case of Jacobi operators H, assuming Hypothesis 2.1. 

Insertion of the norm convergent ansatz 

oo 

m±(z,k) = m± t j(k)z~ j , m± t i{z, k) = -I m (3.1) 
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into (2.41) and (2.42) then yields the following recursion relation for the coefficients 
m± tj (k): 



m+,i = -I m , m+,2 = -B, 
j'-i 

m +d+i = Bm + .j - ^2 Am+ j_ e Am + ,i, j > 2. 



(3.2) 



t=i 



and 



(3.3) 



m_ i = -I m , m_ 2 = ~B, 
j-i 

m-,j+i = Bm-.j — A~mZ j_ e A~rri-,i, j > 2. 
l=i 

Next, rewriting (2.17) in the form 

A{ky 1 M ± {z,k+l)A{k)- 1 M ± {z 1 k) (3.4) 
+ A(k)~ 1 (zl m - B(k + l))A(k)- 1 M ± (z, jfe) + I m = 0, z g C\K, k e Z, 
and inserting the norm convergent ansatz 

oo 

M + (z,k) = J2 M +^( k >~ j (3-5) 
J=l 

and the asymptotic ansatz 

oo 

M_(z,k) = -I m z + *y M- {k)z-' J (3.6) 

into (3.4) then yields the following recursion relation for the coefficients M±j(k): 
M +) i = -A 2 , M +j2 = -AB+A, 

J- 1 (3 7) 

M+j+i = AB+A^M+j - AM+ j _ e A- 1 M +! e, j>2 

and 

M_,o = B, M_,i = (A - ) 2 , M_ i2 = A~B~A~, 

M_ J+1 = -B(A~)~ 1 MZjA~ + ^2M_j_ e (A~)~ 1 MZ ii A~, j > 2. 

e=o 

Remark 3.1. In the continuous cases of Schrodinger and Dirac-type operators dis- 
cussed in detail in [10] and [11], establishing the existence of appropriate asymptotic 
expansions was a highly nontrivial endeavor. Here in the discrete context, hypoth- 
esis (2.2) together with (2.30) immediately yields a norm convergent expansion of 
m±(z, k) as \z\ — > oo. By (2.39) and (2.40), this immediately yields the existence of 
an asymptotic expansion for M±(z, k) as \z\ — > oo. By inspection, these expansions 
are of the form (3.1), (3.5), and (3.6). 

Given the asymptotic expansions (3.5) and (3.6) for M±(z, k) as \z\ — > oo, one can 
of course derive analogous asymptotic expansions for the 2m x 2m Weyl-Titchmarsh 



TRACE FORMULAS AND BORG-TYPE THEOREMS 



11 



matrix M(z, k) in (2.43)-(2.47). For the (1, l)-block matrix element of M(z, k) one 
obtains the norm convergent expansion for \z\ sufficiently large, 



g(z,k) = G(z,k,k) = M 11 (z,k) = V Tj{k)z~ j 

- -Imz- 1 - B(k)z- 2 - [A(k - if + A{kf + B(k) 2 ]z- 3 

- [B(kf + A(k - l)B(k - l)A(k - 1) + A(k)B(k + l)A(k) 



(3.9) 



|z|^OC 



+ B(k)A(k) 2 + B(k)A(k - l) 2 + A(k) 2 B(k) + A(k - l) 2 B(k)] z - 4 
+ 0( Z - 5 ), fceZ. (3.10) 



Similarly, 



G(z,k,k+ 1) 
G(z,k + l,k) 



-A(k) z - 2 + 0(\z\- 3 ), keZ, 
-A(k) z - 2 + Q(\z\- 3 ), keZ. 



(3.11) 
(3.12) 



Next, we also recall that the (1, 1) and (2, 2)-block matrices of M(z, ko) are mxm 
Herglotz matrices. In particular, in addition to (2.48), (2.49), (2.52), and (2.53) 
one obtains 



g(z,k) = G(z,k,k) = M hl (z,k) 

= K 1A (k)+ [ dn 1A (\,k) 

d\E(X,k) 



1 



A 



exp 



1 



A - z 1 + A 2 
A 



z e 



where 
and 



A - z 1 + A 2 



< S( • , k) < I m a.e. 



S(A, k) = lim -Im(ln(.g(A + is, k))) for a.e. A £ R. 

ej.0 7T 



k e Z, 



(3.13) 
(3.14) 

(3.15) 
(3.16) 



We note that the constant term in the exponent of (3.14) vanishes because of the 
asymptotics (3.10). 



4. Borg-Type Theorems for Matrix- Valued Jacobi Operators 

In this section we prove trace formulas and the discrete analog of Borg's unique- 
ness theorem for matrix- valued Jacobi operators. 

In the following, <r(T) and <r css (T) denote the spectrum and essential spectrum 
of a densely defined closed operator T in a complex separable Hilbert space. 

We start with the case of matrix-valued self-adjoint Jacobi operators H assuming 
the following hypothesis: 

Hypothesis 4.1. In addition to Hypothesis 2.1 suppose that o-(H) C [E_,E + ] for 
some — oo < E_ < E + < oo. 

Assuming Hypothesis 4.1, trace formulas associated with Jacobi operators then 
can be derived as follows. First we note that (3.9) implies the expansion (convergent 
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for \z\ sufficiently large) 



d 



- ln(g(z,k)) = fceZ, 
az | z |^oc z — ' 

Si(fc) = J m , 
s 2 (fc) = 

s 3 (fc) - 2A(k - If + 2A(k) 2 + B(k) 2 , etc. 
Moreover, by (3.14), 

d f E+ f°° 

—ln(g(z,k)) = / d\(\- z)- 2 ~(\,k)+ / dX(X-z)- 2 I„ 
dz Je_ J e + 

= (E+ - z)-H m + d\ (A - z)- 2 S(A, fc), 



(4.1) 

(4.2) 
(4.3) 
(4.4) 



where we used 



0, A<£_, 
I m , A > 



H(A,fc) = 

Theorem 4.2. Assume Hypothesis 4-1- Then (cf. (4.1)), 



(4.5) 



(4.6) 



i(k) = ±(EL- 1 + E*r 1 )I 7n + ±(j-l)J + d\\i- 2 [I m -2Z(\,k)], j en, fcez. 



Explicitly, for all fceZ, 
S(A ; ) = i( J B_+ J E; + )7 ro + i 



(4.7) 
(4.8) 



dX [I m -2E(X,k)}, 

2A{k-l) 2 + 2A{k) 2 +B{k) 2 = \{E 2 _ +E 2 + )I m + j + dX X[I m - 2 H(A, fc)], etc. 

(4.9) 

Proof. By (4.1) and (4.5) one infers 



z-E+ 



E. 



|z| — »OC 



dX (X-z)- 2 [I m -2~(X,k)} 
(4.10) 
(4.11) 



Expanding (4.10) with respect to z 1 and comparing powers of z J then yields 
(4.7). Equations (4.8) and (4.9) arc then clear from (4.3) and (4.4). □ 

Remark 4.3. In the scalar case m = 1, Theorem 4.2 was first derived in [26] 
assuming A(k) = 1, k £ Z. The Jacobi case for half-lines was explicitly discussed 
in this vein in [28]. The trace formula (4.8) for full-line Jacobi operators in the case 
m = 1 (and other trace formulas) can be found in [53], [54, Sect. 6.2]. The current 
matrix- valued trace formula for to > 2 is new. 
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Next we turn to a Borg-typc theorem for matrix- valued Jacobi operators. To 
set the stage we first recall Flaschka's result [19] for scalar (i.e., m = 1) periodic 
Jacobi operators, the direct analog of Borg's theorem for periodic one-dimensional 
Schrodinger operators originally proved in [7]. 

Theorem 4.4. ([19]) Suppose a and b are periodic real-valued sequences in £°°(Z) 
with the same period and a(k) > 0, k E Z. Let h(a, b) = aS + + a~ S~ + b be the 
associated self-adjoint Jacobi operator in l 2 (1) (cf. (2.4) for m = 1) and suppose 
that a{h{a, b)) = [E_, E + ] for some E_ < E + . Then, 

a{k) = {E + -E_)/A, b(k) = {E- + E+)/2, fceZ. (4.12) 

While uniqueness results such as Theorem 4.4 are described as Borg-type theo- 
rems, other types of inverse spectral results are also described as Borg-type theorems 
as mentioned for Schrodinger operators in Remark 1.2. We also note that Theorem 
4.4 is quite different from a recent result by Killip and Simon [34, Thm. 10.1] which 
states that a(k) = 1, k E Z and a(h) C [—2, 2] implies b(k) — 0, k E Z. 

As shown in [11] and [13], periodicity is not the key ingredient in Borg-type 
theorems such as Theorem 4.4. In fact, it was shown there that the more general 
notion of being rcflcctionless is sufficient for Borg-type theorems to hold and we 
will turn to this circle of ideas next. We note that the class of rcflcctionless in- 
teractions include periodic and certain cases of quasi-periodic and almost-periodic 
interactions. 

Following [11] and [13] , we now define reflectionless matrix- valued Jacobi oper- 
ators as follows: 

Definition 4.5. Assume Hypothesis 4.1. Then the matrix-valued coefficients A, B 
are called reflectionless if for all k £ Z, 

S(A, fc) = for a.c. A E a ess (H) (4.13) 

with S(-,fc) defined in (3.16). 

Since hardly any confusion can arise, we will also call H = H(A, B) reflectionless 
if (4.13) is satisfied. 

Remark 4.6. In the next theorem we will prove an inverse spectral result for 
matrix- valued Jacobi operators H(A,B). However, in the general matrix- valued 
context, where m > 2, one cannot expect that the spectrum of H (A, B) will de- 
termine A and B uniquely. Indeed, assume that B is a multiple of the identity, 
B(k) = b(k)I m for some b E r°(Z), b(k) E R, k E Z. In addition, let U be a unitary 
mxm matrix and consider A(k) = UA(k)U~ 1 , k E Z. Then clearly H(A,B) and 
H(A, B) are unitarily equivalent and hence the spectrum of H cannot uniquely 
determine its coefficients. The following result, however, will illustrate a special 
case where the spectrum of H(A, B) does determine A and B uniquely. 

Given Definition 4.5, we now turn to a Borg-type uniqueness theorem for H 
and formulate the analog of Theorem 4.4 for reflectionless matrix-valued Jacobi 
operators. 

Theorem 4.7. Assume Hypotheses 4-1 and suppose that A and B are reflectionless. 
Let H(A, B) = AS + + A~S~ + B be the associated self-adjoint Jacobi operator in 
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£ 2 {Z) m (cf. (2.4)) and suppose that a(H(A, B)) = [E-,E+] for some E- < E+. 
Then A and B are of the form, 



A(k) = -{E+ - E_)I n 



B(k) = -(£_ + E+)I„ 



fc g z. 



(4.14) 



Proof. By hypothesis, S(A, fc) = (l/2)7 m for a.e. A G [E_,E+] and all k G Z. Thus 
the trace formula (4.8) immediately yields (4.14) for B. Inserting the formula (4.14) 
for B into the second trace formula (4.9) one infers 

A(k-1) 2 + A(k) 2 = l(E+-E_) 2 I m , keZ. (4.15) 



The first order difference equation (4.15) has the solution 



A(2£) 2 = A(0) 2 , A(2£ + l) 2 = \{E+ - E_) 2 I n 

8 



A(0) 2 , £eZ. (4.16) 



Since by hypothesis A(0) > is a self-adjoint m x m matrix, there exists a unitary 
mxm matrix U that diagonalizes A(0) and by (4.16), U simultaneously diagonalizcs 
A(k) for all keZ, 

A(k) = UA(k)U~ 1 , keZ, (4.17) 

where A(k) are diagonal matrices for all k e Z. By (3.5)-(3.8) followed by an 
analytic continuation to all of C+, U also diagonalizes M±(z, fc), 

M±{z,k) = UM±(z,k)U- 1 , fceZ, (4.18) 

where M±(z, fc) are diagonal matrices for all fc G Z. The same result of course fol- 
lows from (2.16) taking into account that U also diagonalizcs ip±i z i k, ko), 8(z, k, fco), 
and 4>(z,k,ko)- The resulting diagonal matrices will of course be denoted by 
ip±(z, k,ko), 9(z,k,ko), and 4>(z,k,ko) below. 

Next, we will invoke some Herglotz function ideas. Since for all fc G Z, 



0, 



5(A,fc) 



2 In 



A < E-, 
A G (E_ , E + ) 

Im, A > E + , 

we can compute g(z 7 k) from (3.14) and obtain 

il\ - ( A. A ) | — 

*E 



(4.19) 



g(z, fc) = exp 



exp 



r E + 
(1/2) / 

+ / dX 



+ \ 2 

A 



A - z 1 + A 2 
A 



X-z 1 + A 2 , 

= -C[(z - E-){z - E+)]- l/2 I m , zgC+, fcGZ (4.20) 

for some C > 0. However, the known asymptotic behavior (3.10) of g(z,k) as 
\z\ — > exo then proves C = 1 and hence 

g(z,k) = -{(z-E_)(z-E + )}- 1 / 2 I m :=g(z), z G C+, fc G Z. (4.21) 

Next, we consider —g(z)~ 1 and its Herglotz representation, 

-giz)- 1 - M+(*, fc) - M_(z, fc) = [(z - E_)(z - E + )] x / 2 I m 
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= zI m -B + [ dT(X)(X-z)- 1 

JR 

= f dfi+(A,fc)(A- z)- 1 + zl m - B+ [ dQ_{X,k)(X- z)' 1 , (4.22) 

JR JR 

where dCl±(X, k) denote the measures in the Herglotz representations of ±M±(z, k), 

M + (z,k)= [ dn + (\,k) (A - z)- 1 , (4.23) 
Jr 

-M-(z,k) = zI m -B+ [ dn_(A,jfe)(A-z) _1 . (4.24) 

Jr 

Actually, applying U from the left and U^ 1 from the right on either side in (4.22)- 
(4.24), we may replace M± and dT, dfl± by M± and dT, dQ.± (where in obvious 
notation dT and dCt± denote the diagonal measures associated with the diagonal 
Herglotz matrices M + — M_ and M±, respectively). In the following we assume 
that all quantities in (4.22)-(4.24) have been replaced by their diagonal matrix 
counterparts. 

Next, we follow a strategy employed in [24] in the scalar Jacobi case. First we 
note that S(A, k) = (1/2)7 TO for a.c. A E (E_,E + ) is equivalent to 

g(\ + iO) = -g(X + iO)*, A e (£_,£+) (4.25) 

and hence to 

-giX + iQ)- 1 = [g(\ + iQ)- 1 }*, A £(£_,£+). (4.26) 

Here and in the following, /(A + iO) denotes the normal limit lim e | /(A + ie). The 
last result is easily seen to be equivalent to the following fact: For all k E Z, 

Re(M+(A + iO,fc)) = Re(M_(X + iO,k)) for a.c. A E (E_,E + ). (4.27) 

By (2.47) and (2.55) one computes 

UA(k)g(z, k + l)^)!/" 1 = A(k)g(z, k + l)A(k) 

= UM±(z,k)[M-(z,k) - M + (z,k)]- 1 M T (z,k)U- 1 

= M±(z, k)[M-(z, k) - M+(z, k^M^z, k) 

= M±(z,k)g(z,k)M T (z,k), zeC+, fceZ. (4.28) 

Since g(z, k) — g(z) is independent of k E Z, (4.28) implies 

A(k) 2 =M + (z,k)M^(z,k) =M-(z,k)M+(z,k), z E C+, k E Z. (4.29) 

Inserting M±(z,k) = Re(M±(z,k)) + ilm(M±(z, k)) into (4.29) then explicitly 
yields 

lm(M + (z,k))Re(M-(z,k))+Re(M + (z,k))Im(M_(z,k)) = (4.30) 

and since M±(z, k), Re(M±(z, k)), and lm(M±(z, k)) are all diagonal matrices, we 
note that all entries in (4.30) commute. Similarly, the fact 

g(X + iO,k + l) 

= A{k )- 1 M±{X + iQ 1 kv) 

x [M-(X + iO,k a ) -M + (A + iO,fc )] _1 M T (A + iO,A:o)l(fco) _1 
= -g(X + iO, k Q + 1)* = -i(fc ) _1 M T (A + *0, fc )* 
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x [M-(\ + iO,k )* -M + (A + iO,/co)*] _1 M ± (A + iO,/c )*I(fco) _1 (4.31) 

for a.c. A G (E-,E+) 

implies 

M±(A + iO, fc )[M_(A + iO, fc ) - M+(A + iO, fc )]" 1 M T (A + iO, fc ) 

= M T (\ + iO,k a )*[M-(\ + iO,k a ) -M + (A + iO,fc )] _1 M ± (A + iO,/co)* (4.32) 

for a.c. A G (E_,E + ) 

since by (4.21), 

5 (z,fc) = [M.^^-M+^fco)]" 1 - [M_ (*,*<>) -M+(z,*b)] _1 (4-33) 

and hence (4.25) also applies to [M_(A + iO, fc ) - M + (A + iO, fco)] -1 , A G (£_,£+). 

Inserting expression (2.8) for -0± in terms of 0, 0, and M± into (2.38) taking 
£ = k, and inserting the result into (4.26) then yields, 

- (f>(X, fc, k )A{k )- 1 M±(\ + iO, fc )[M_(A + iO, fc ) - M + (A + iO, fco)]- 1 

x 6>(A,fc,fc )* 

- 6(\, k, fc )[M_(A + iO, fco) - M+(A + iO, fc )]" 1 M T (A + iO, fco)A(fco)" 1 

x </>(A, fc, fco)* 

= -<£(A, fc, fc )A(fc )- 1 M T (A + iO, fc )* [M_(A + iO, fc ) - M+(A + *0, fc )] _1 
x 6»(A,fc,fc )* 

- <9(A, fc, fc )[M_(A + iO, fc ) - M+(A + iO, fc )] _1 M ± (A + iO, fco)*A(fc ) _1 
x0(A,fc,fc o )*, (4.34) 

where we also used (4.25) and (4.32). Applying U and U~ x from the left and right 
on either side in (4.34), using the fact that 

4>{z,k,k )* = 0(z,fc,fc o ), 9(z,k,k )* =6(z,k,k ), (4.35) 

and that all diagonal matrices commute, one can rewrite (4.34) in the form, 

2i4>{\, fc, fc o )0(A, fc, fco)A(fco) _1 5(A + iO) (4.36) 

x [Im(M_(A + iO,fc )) + Im(M + (A + iO,fc ))] = for a.e. A G (E_,E+). 

Since fc G Z can be chosen arbitarily in (4.36), this implies 

Im(M_(A + iO,fc )) = -Im(M+(A + iO,fc )) for a.e. A G (E-,E + ). (4.37) 

Finally, since fco G Z in (4.37) is arbitrary, one obtains for all fc G Z, 

Im(M_(A + iO,fc)) = -Im(M+(A + iO,fc)) for a.e. A G (E_,E + ) (4.38) 

and hence also for all fc G Z, 

Im(M_(A + iO, fc)) = -Im(M+(A + iO, fc)) for a.e. A G (E_,E + ) (4.39) 

applying I/ -1 and U from the left and right on both sides in (4.38). Together with 
(4.27) this yields for all fc G Z, 

M_(A=F«0,fc) = M_(A±iO,fc)* = M+(A±iO,fc) for a.e. A G (E-,E + ). (4.40) 
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Thus, M_(-,fc) is the analytic continuation of M + (-,k) (and vice versa) through 
the interval (E-,E + ). Since dT is purely absolutely continuous (cf. (4.22)), 

dT = dT ac , cffpp = rfr sc = (4.41) 

(where d[i ac , d/j pp , and dfj, sc denote the absolutely continuous, pure point, and 
singularly continuous parts of a measure d/i), one also infers 

dQ ±;PP = dfi±, sc = 0. (4.42) 

(This also follows from the fact that M±(-,k) have analytic continuations through 
(E-,E + ), see [30, Lemma 5.6].) Especially, (4.38) then implies the ^-independence 
of dQ±(-, k), 

dCl+(X, k) = dft_(A, k) = ^dT(A). (4.43) 
Equations (3.7), (3.8), (4.23), (4.24), and (4.43) then prove 



A(k) 2 = f + dn+(X,k) = \ f + dT(\)= f + dfi_(A,fc) = A{k-l) 2 , fceZ 
JE- 2 J E _ J E _ 

(4.44) 

and hence also 

A(k) 2 = A{k- l) 2 , fceZ, (4.45) 
which proves that A is independent of k G Z. The trace formula (4.15) for A, 

A(k-l) 2 + A(k) 2 = \{E + -E_) 2 I m , keZ. (4.46) 

o 

then proves 

A(k) 2 = -^(E + - E_) 2 I m , fc e Z. (4.47) 
Using (4.17) then completes the proof of (4.14). □ 
Because of (4.43) and 

dm = S^- E -^~^ lelE-.*], (4 . 48) 
[0, A e M\[E-,E+], 

(cf. (4.22)), one obtains 

M ± (z,k) = ^-\z + \{E_ + E + )±^[{z-E_){z-E + )] 1 ^I m , zeC+ , k e Z. 

(4.49) 

Corollary 4.8. ^4sswme Hypothesis 4-1 in the special case m = 1 and suppose that 
a and b are reflectionless. Let h(a,b) — aS + + a~ S~ + b be the associated self- 
adjoint Jacobi operator in £ 2 (Z) (cf. (2.4)) and suppose that a(h(a,b)) = [E-,E + ] 
for some E_ < E + . Then, 

a(k) = ^(E + -E_), b(k) = ^(E.+E+), fceZ. (4.50) 

While Theorem 4.7 is new, Corollary 4.8 in the scalar case m = 1 was noted in 
[53, Corollary 6.3] (see also [54, Corollary 8.6]). 

The following result can be proved in analogy to Theorems 4.6 and 4.8 in [13], 
hence we state it here without proof. 
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Theorem 4.9. In addition to Hypothesis J^.l, suppose that A and B are periodic 
with the same period. Let H(A,B) = AS + + A~S~ + B be the associated self- 
adjoint Jacobi operator in £ 2 (Z) m and suppose that H(A,B) has uniform spectral 
multiplicity 2m. Then H{A, B) is reflectionless and hence for all k e Z, 



In particular, assume that A and B are periodic with the same period, that H( A, B) 
has uniform spectral multiplicity 2m, and that a(H(A, B)) = [E_,E + ] for some 
E_ < E + . Then A and B are of the form, 



In connection with the special case m = 1 in Theorem 4.4 we note that scalar Ja- 
cobi operators automatically have uniform spectral multiplicity 2 since the product 
of the two Floquet multipliers equals one. 



In our final section we turn to a Borg-type theorem for a class of supersymmetric 
Dirac difference operators. Rather than developing the results from first principles 
as in the case of Jacobi operators with matrix-valued coefficients, we will employ 
the underlying supersymmetric structure and reduce the case of Dirac difference 
operators to that of Jacobi operators. 

We start with the following abstract result. 

Theorem 5.1 ([14], [25]). Let C be a densely defined closed operator in a separable 
complex Hilbert space H with domain V(C) and introduce the operator 



S(A,fc) = - 



I m for a.e. A E cr ess (H(A, B)). 



(4.51) 



A(k) = \{E+- E_)I m , B(k) = \{E-+ E+)I m , k e Z. 



(4.52) 



5. Borg-Type Theorems for Supersymmetric Dirac 
Difference Operators 




(5.1) 



inH®H. Then, 

Q = Q*, 



(5.2) 




(5.3) 



(5.4) 



(5.5) 



z 2 e C\{a(C*C)Ucr(CC*)}. 




(IhT 1 2 C(C*C - c/w) _1 c*, 

C G C\{a(C*C) Ua(CC*)}, 

c^r 1 2 c*(cc* - c/w) _1 c, 

C eC\{a(C*C)Ua(CC*)}. 



(5.6) 



In + C(C*C 



(5.7) 
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Moreover, 

implies 
and hence 
Conversely, 
implies 

and 
implies 



QU(z) = zU(z), U(z)=(^j (5.8) 

C*u 2 (z) = zu\(z), Cu\(z) — zu 2 (z) (5.9) 
C*C Ul (z) = z 2 Ul (z), CC*u 2 (z) = z 2 u 2 (z). (5.10) 
C*Cui{z) = z 2 Ul (z), z^O, (5.11) 

QU(z) = zU(z), W=( (1/ ^ w ) (5-12) 
CC*u 2 (z) = z 2 u 2 {z), z^O, (5.13) 

QU(z) = zU(z), U(z) = ( {1/Z) U C 2 ^ 2{Z) ) • (5-14) 

In order to apply this setup to finite difference Dirac-type systems (cf. [12]), we 
introduce the following hypothesis. 

Hypothesis 5.2. Let m e N and consider the sequence of invertible mxm matrices 
P = {p(k)}kez e C(Z) mxm , p(k) = p(k)*, keZ, (5.15) 
X = {xW}«eC(Zf xm , (5.16) 
det C m(p(fc)) ^0, det C m(x(fc)) ^ 0, fcel (5.17) 

In addition, we assume that p{k) is a positive definite diagonal mxm matrix 

/ o(fc) = diag(pi(fc),...,p m (fc)), pj(k) > 0, 1 <j < m, k e Z, (5.18) 

and that p(k) and %(fc) are uniformly bounded with respect to k £Z, t/iat is, i/iere 
exists a C > 0, suc/i i/iai 

||p(*0llc™*™ + llxWIIc™*™ <C, fcGZ. (5.19) 

Finally, we suppose that p\ + and \p are positive definite, 

p(k) X (k + l) >0, xWp(fc) > 0, fc e Z. (5.20) 

Assuming Hypothesis 5.2, we thus introduce the bounded linear operator 

E = p-S-+ X , V(E)=£ 2 (Z) m , (5.21) 
on £ 2 (Is) m and the bounded Dirac-type difference operator 

D= (e E ()) =S p + X > V(D)=l 2 (Z) rn (Bl 2 (Z) m (5.22) 
on £ 2 (Z) m ®£ 2 (Z) m , where 



o P s+\ (0 x 



s ?={ps- o J' x =[x o)- ^ 

One then computes 

H 1 = E*E = A 1 S + +A^S~+B U (5.24) 
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H 2 = EE* = A 2 S+ + A^S- +B 2l (5.25) 

where 

A 1 (k) = p(k) x + (k) > 0, B 1 (k) = ( P (k)f + X (kYx(k), keZ, (5.26) 

A 2 (k) = X (k)p(k) > 0, B 2 (k) = (p(k)-) 2 + x(k) X (kr (5.27) 

and notes that Hi > and H 2 > are matrix- valued Jacobi operators in £ 2 (Z) m 
of the form (2.4). 

Remark 5.3. We note that D with a positive definite diagonal m x m matrix p 
in (5.22) represents a normal form of Dirac-type difference in the following sense: 
Assume Hypothesis 5.2 except for the condition that p(k) is a positive definite 
diagonal matrix for all k G Z. Then, following [12, Lemma 2.3], there exists a 
sequence of unitary matrices U(p) = {U(p,k)}kez G C(Z) 2mx2m such that 

U(p)(S p + X)U(p)- 1 =Sp + X, (5.28) 

where p is diagonal and positive definite and X is of the form 

X = U(p)XU( P y 1 = Q (5.29) 

with x G C(Z) mxm . Indeed, denote by Q p {k) G C mxm a unitary matrix such that 
Qp(k)p(k)Q p (k)^ 1 = p(k), where G R mxm is diagonal and self-adjoint for all 
keZ. Then, 



(7 P (S P + X)^ 1 = S~ p + X, X = UpXU- 1 , Up = (Of Q °J . 



(5.30) 



Next, let e{k) G R mxm be a diagonal matrix for which {e(k)) e ,e G {+1, -1}, I = 
1, ... , m. Define e(fc) G K mxm by e(fc) = e(k)e(k + 1) and choose e(fc) so that 
p = ep > 0. Then, 

£40% + A^C/- 1 = ,Sp + X, X = U £ XU-\ U £ = °) . (5.31) 
Thus, one obtains 

U(p) = U s U p =^ ^_), (5.32) 

X = Q , x = eQ'xQ^e. (5.33) 

Next, we introduce the 2m x 2m Weyl-Titchmarsh matrices associated with D 
(cf. [12]) by 

M D (z,fc) = (M^,(z,fc)) w . /=12 , zeC\R, fceZ, (5.34) 

where 

M[> tl (z,k) = [M D (z,k) -M°(z,k)]~\ (5.35) 

M® 2 (z, k) = 2- 1 [M D (z, k) - M°(z, k)] _1 [M D (z, k) + M°(z, k)] , (5.36) 

M^(z, k) = 2- 1 [M D (z, k) + M°(z, k)] [M D {z, k) - M»(z, k)] ~\ (5.37) 

M° 2 (z,k) = M£(z,k)[M D (z,k)-M°{z,k)]~ 1 M T (z,k), (5.38) 
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and similarly, the 2m x 2m Weyl-Titchmarsh matrices associated with the Jacobi 
operators H e , £ = 1,2 (cf. (2.43)-(2.47)) by 

M H <(z,k)=(M*f,(z,k)) jt . l=U2 , *=1,2, *eC\R, fceZ, (5.39) 

where 

M*{{z,k) = [M He (z,k)- M^(z,k)]~\ (5.40) 

M$ (z, k) = 2- 1 [M He (z, jfc) - Ml ' (z, jfe)] _1 [M*< (z, jfc) + (z, jfe)] , (5.41) 

M*[ (z, jfc) = 2- 1 [M He (z, jfe) + Aff f (z, jfe)] [M Hf (z, jfe) - (z, jfc)] ~\ (5.42) 

M%(z,k) =Mg t (z,k)[M? t (z,k)-Ml l (z,k)]~ 1 M T (z,k). (5.43) 

The supersymmetric formalism (5.1)-(5.5) then implies the following relations be- 
tween M D (z) and M He (z), £=1,2. 

Theorem 5.4. Assume Hypothesis 5.2 and let z € C\a(D), k G Z. Then, 

Mg(z,k) = -z-V(fc) + ^ _ V(fc) _1/2 ^f 1 (2 2 ,fc)p(fc)- 1 / 2 , (5.44) 
Mj?(z, fc) - -zp(k)- 1 - zpik)- 1 ' 2 [xiky'M^ixik)- 1 )*} ^p(k)- 1 / 2 . (5.45) 

Proof. We freely employ the relations (5.8)-(5.14) and some results from [12]. Let 

U±(z,k,ko)=h^l>t°l), u l!± (z,k ,k )=I m (5.46) 



K U2,±(z, k, A*)), 

be the normalized Weyl-Titchmarsh solutions associated with D. Then by equa- 
tions (2.35a)-(2.35c) and (2.96) in [12], 

u 2 ,±{z, k + 1, k )u h ±{z, k, ko)- 1 = -p(k)-^ 2 M^(z, k)p(k) 1 / 2 . (5.47) 

Moreover, using (C*U2 : ±)(z, ko, k ) — zu\.±{z, k , k ) = zl m (cf. (5.9)), one derives 

u 2 . ± (z, k , k ) = (xikoTr 1 [zl m + p(k ) 1/2 M°(z, k )p(k ) 1/2 ] . (5.48) 
Next, let 

tpe,±(z,k,ko), ipe,±(z,ko,ko) = I m (5.49) 
be the normalized Weyl-Titchmarsh solutions associated with H e , I = 1,2. Then 
by (2.16), 

Mi e (z,k) = -A t {k)^ ± {z,k+l,ko)i)i,±(z,k,ko)- 1 . (5-50) 

Given the uniqueness of Weyl-Titchmarsh solutions for D and Hi, I = 1,2, (5.8)- 
(5.14) yield 

ui,±(z, k, k Q ) = ip li± (z 2 ,k, ko), (5.51) 

u 2 ,±(z,k,ko) = (l/z)Eui i± (z,k,k ) 

= (l/z)[p~(k)ui t ±(z, k-1, k ) + x( k )ui,±(z, k, k )] 

= (l/z)[p-(k)^ h± (z 2 ,k - 1, ko) + X(k)ipi,±(z 2 ,k, k )]. (5.52) 

Thus, 

u 2 ,±(z, k+1, ko)u h ±{z, k, ko)- 1 = -p{k)-^ 2 M^{z, k)p(k)^ 2 

= (1/z) [p{k)^±{z 2 , k, k ) + X (k + l)ipi,±(z 2 , k + 1, ko)] <K±(z 2 , k, ko)- 1 
= (1/z) [p(k) + X (k + lH-A^k))-^ (z 2 , k)} 
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= (l/z)[p(k)-p(k)- 1 M^(z 2 ,k)] (5.53) 
proving (5.44). 

Similarly, the uniqueness of Weyl-Titchmarsh solutions also yields 

u 2 ,±(z, k, fc ) = i/> 2 ,±(z 2 , fe, fco)d±(z, fe ) (5.54) 
for some constant to x to matrix d±(2, fco). Thus, 

d±(*,fco) =«2,±(«,fco,*o) - (#o)T 1 K + #o) 1/ X («,*o)p(*d) 1/2 ]- (5.55) 
One then computes, 

U2,±(«, fco + 1, fco) - -p(fc )- 1/2 M| ) (z, k )p(k ) 1/2 

= ^ 2 ,±(^ 2 , fc + 1, fc )(x(fc )T V™ + p(Ao) 1/2 Mf («, fc )p(fc ) 1/2 ] 

= -(A 2 (fco))- 1 Mf(z 2 ,fco)(x(fco)T 1 [^ + #o) 1/2 M ± D (z,io)#o) 1/2 ] 

= -p(fco)- 1 x(fco)" 1 Mf 2 (z 2 ,fc )(x(fco)*)" 1 [^m + p(fco) 1/2 M| ) (z,fco)p(fco) 1/2 ]. 

(5.56) 

Hence, 

p(fc ) 1 /2 M | , (z,fc )p(A :o ) 1 /2 

= x(fco)- 1 Mf( z 2 ,t )( x (fc )T 1 ^ m +p(^) 1/2 M ± D (^o)p(fco) 1/2 ] (5.57) 
and since ko £ Z is arbitrary, 

p{k) 1 ' 2 M^{z,k)p{k) 1 ' 2 (5.58) 
= X (fc)- 1 Mf 2 (z 2 ,A : )(x(A : )*)- 1 [^ m + pW 1/2 M| , (z,A : )p(fc) 1 /2] j fceZ . 

Solving (5.58) for Af£(z, fc) then yields (5.45). □ 

According to (2.49), M D (z 1 k) is a matrix-valued Hcrglotz function of rank 2to 
with exponential Hcrglotz representation 



M D (z 1 k)=cxp\C' D (k) + J^d\? D (\,k)\j— - T -^jj, (5.59) 

where 

C D (fc) = C D (£;)*, 0«:T D (.,fc)^/ 2m a.e, (5.60) 

T D (\, k) = lim -Im(ln(M(A + is, k))) for a.e. A e R. (5.61) 

Following [11] we now define reflectionless matrix- valued Dirac-type operators 
as follows: 

Definition 5.5. Assume Hypothesis 5.2. Then the matrix-valued coefficients p, x 
are called reflectionless if for all k £ Z, 

T D (A, fc) = i/ 2m for a.e. A e a ess (D). (5.62) 

We also call D = D(p,\) reflectionless if (5.62) holds. 

Remark 5.6. The definition (4.13) of reflectionless Jacobi operators H and the 
definition (5.62) of reflectionless Dirac operators D can be replaced by the more 
stringent requirement that for all k £ Z, 

M+(A + iO, k) = M_(A - iO, k) for a.e. A £ a css (H) (5.63) 
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for Jacobi operators H and similarly for Dirac-type operators D replacing M±(z, k) 
by M± (z, k). This yields a unified definition of the notion of rcflcctionless matrix- 
valued Jacobi and Dirac-type operators. It is easy to see that (5.63) implies (4.13). 
The converse is more subtle and was proved by Sodin and Yuditskii in [51] and [52] 
for scalar Schrodinger and Jacobi operators H under the assumption that a(H) is 
a homogeneous set. Their proof extends to the present matrix- valued setting. 

In the special case of a Borg-type situation with cr(H) = [E-,E + ], we explicitly 
derived (5.63) in (4.40). In this particular case, M_(-,fc) is the analytic continu- 
ation of M + (-,k) (and vice versa) through the interval (E_,E + ). In general, the 
homogeneous set a(H) may be a Cantor set (of positive Lebesguc measure) and 
then M±(-, k) are called pseudo-continuable through <j(H). 

Lemma 5.7. Assume Hypothesis 5.2. If D is reflectionless (in the sense of (5.62)) 
then Hi, i = 1,2, are reflectionless (in the sense of (4.13)). 

Proof. D being reflectionless in the sense of (5.62) is equivalent to the assertion 
that for all k G Z, M D (X + iO, k) is skew-adjoint for a.c. A e a css (D). Equivalently, 
for all k G Z, 

M D (X + iO, k) = iC(X, k) with C(A, jfe) = C(X, k)* for a.e. A e cr css (L>). (5.64) 

In fact, C(A, k) — lm(M D (X + iO, k)) > 0. As a consequence, also all block subma- 
trices of M D (X + iO, k), symmetric with respect to the diagonal of M D (X + iO, k), 
are skew-adjoint. In particular, the two m x m diagonal blocks of M D (X + iO, k) 
satisfy for all k e Z, 

Mfj(X + iO, k) = iC\ e (X, k) with Co(A, k) > for a.e. A e a ess (D). (5.65) 

Equation (5.44) implies 

[M D (z, k) - M°(z, k)}- 1 = p(k) 1 / 2 z[M Hl (z 2 , k) - Ml 1 (z 2 , k)] '^(k) 1 / 2 , fceZ 

(5.66) 

and hence, for all k G Z and a.e. A G a css (D), 

M^(X + iO, k) = [M D (X + iO, k) - Mf (A + iO, k)] ^ 

= p(k) 1 / 2 X[M Hl (A 2 + iO, k) - Ml 1 (A 2 + iO, k)\ '^(k) 1 / 2 . (5.67) 

Thus, Mfl (X 2 +i0, k) = g Hl (X 2 +i0, k) is skew-adjoint and hence Hi is reflectionless 
by (4.13).' 

Similarly, applying (5.45) yields 

z[M H2 (z 2 ,k)-Ml 2 (z 2 ,k)]~ 1 

= (x(k)- 1 Y p(k) 1 ' 2 [M D (z, k) + zp(k)- 1 ] [M D (z, k) - M?(z, k)] - 1 

x [M?(z,k) + zp(k)- 1 ]p(k) 1 / 2 x(ky 1 , fceZ (5.68) 

and hence, for all k £ Z and a.e. A G a css (D), 

X [M H2 (A 2 + iO, k) - Ml 2 (A 2 + iO, k)] ^ 

= ( X (k)- 1 T P (k) 1 ' 2 [M D (X + iO, k) + Xp(k)- 1 ] 

x [M D (X + iO,k)-Ml(X + iO,k)y 1 

x [Ml(X + iO,k) + Xp(k)- 1 ]p(k) 1 / 2 X (k)- 1 . (5.69) 
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Next, consider the Herglotz matrix 
M a {z,k) = [M D (z,k) + a(k)] [M D (z,k)-M 1 j ) (z,k)y 1 [M°(z,k)+a(k)] (5.70) 

for some self-adjoint m x m matrix a(k). We claim that for all k <G Z and a.e. 

A e a ess (D), 

M a (X + iO,k) = iC a (X,k) with C a {X,k) = C a (X,k)*. (5.71) 

Indeed, one computes 

M a (\ + i0,k) = M D (X + iO,k) [M D (\ + iO,k) -M£(A + iO,ft)] _1 -M+ (A + iO, fc) 

+ a(jfe) [M D (A + iO, k) - Mf (A + iO, fc)] _1 a(A:) 

+ a(fc) [M D (A + iO, fc) - M^(A + »0, fc)] _1 m£(A + iO, fc) 

+ Mf(A + iO, fc) [AT? (A + iO, fc) - Mf (A + iO, fc)] _1 a(*0- (5.72) 

The first two terms on the right-hand side of (5.72) are already of the required form 
(5.71). The last two terms on the right-hand side of (5.72) can be rewritten in the 
form 

a(k) [ME (A + iO, k) - Aff (A + iO, fc)] _1 Aff (A + iO, k) 

+ M D (\ + iO, k) [ME (A + iO, fc) - Mf (A + *0, fc)] _1 a(fc) 
= (l/2)a(fc) [M D (A + iO, fc) - M^(A + *0, fc)] _1 

x [ME(X + iO, k) + Mf (A + iO, fc)] 

+ (1/2) [AT? (A + iO, fc) + Aff (A + *0, fc)] 

x [M D (A + iO,fc) -M^(A + iO,/c)] _1 a(/c) 
= (l/2)i[a(jfe)Ci, 2 (A, fc) + Ci, 2 (A, fe)*a(fc)] (5.73) 

using (5.64). Thus, also the last two terms on the right-hand side of (5.72) are 
of the required form (5.71) which completes the proof of the claim (5.71). The 
result (5.71) applied to (5.69) then proves that M^(A 2 + iO, k) = g H2 {X 2 + iO, k) 
is skew-adjoint and hence also H 2 is rcflectionless. □ 

The next result is a Borg-type theorem for supcrsymmctric Dirac difference op- 
erators D. However, unlike the Borg-type theorem for (matrix-valued) Schrodinger, 
Dirac, and Jacobi operators, this analog for Dirac difference operators displays a 
characteristic nonuniqueness feature. 

Theorem 5.8. Assume Hypothesis 5.2 and suppose that p and \ are reflection- 
less. Let D(p,x) — (p-g°- +x P ' S ' +X ) ^ e associated self-adjoint Dirac dif- 
ference operator in £ 2 (Z) m © £ 2 (Z) m (cf. (5.22)) and suppose that a(D(p,x)) = 
[ - E^ 2 , -E 1 ! 2 ] U [e]/ 2 ,E^ 2 ] for some < E_ < E+. Then p and X are of the 
form, 

p(k) = diag(pi(fc), . . . ,p m (k)), 

Pi (k) = \ « /2 - ejEl /2 ) , 1 < j < m, k G Z, (5.74) 
X(fc) = diag(xi(fc), . . . ,Xm(k)), 
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Xj(k) = \ (El /2 + e,£ V2 ) , 1 < 3 < m, k G Z, (5.75) 

where 

ejG{l,-l}, l<j<m. (5.76) 

Proof. Since D is rcflcctionlcss by hypothesis, so are Hi and H 2 by Lemma 5.7. By 
(5-3), 

a(D) = [ E l f, -E 1 ! 2 ] U [E 1 / 2 , E l f] (5.77) 

implies 

a(H e )=[E_,E+], £=1,2. (5.78) 
Applying Theorem 4.7 to Hi and H 2 then yields 

- A {E + -E_)I m = Ai=px + = A 2 = XP , (5.79) 

i(£+ + £L)J m = Bi = p 2 + X *X = B 2 = (p-f + XX*, (5.80) 
By (5.79), x satisfies 

X = X + = \{E+-E_)p- 1 (5.81) 

and hence ^ is a constant (i.e., k- independent) positive definite sequence of m x 
m matrices. By (5.80) (equivalently, by (5.81)), then also p = p~ is a constant 
sequence ofmxm matrices. Insertion of (5.81) into (5.80) then yields 

P 4 - \(E+ + E_)p 2 + ^(E+- E_ fl m = 0. (5.82) 

Since by hypothesis p is a positive definite diagonal matrix, solving the quadratic 
equation (5.82) for pj, 1 < j < m, yields (5.74) and hence also (5.75) using (5.81). 

□ 

To the best of our knowledge this result is new even in the scalar case m = 1. 
In particular, the sign ambiguities displayed in (5.74) and (5.75), giving rise to 2 m 
isospectral supersymmetric Dirac difference operators, are new in this Borg-type 
context. The sign ambiguity disappears and hence uniqueness of the correspond- 
ing inverse spectral problem is restored only in the special case E- = 0, that is, 

/ 1 /2 1 /2\ 

whenever the spectral gap ( — E_ ,E_ ) of D(p,x) closes. 

Using the supersymmetric formalism described in this section, the proof of fol- 
lowing result can be reduced to that of Theorem 4.9. 

Theorem 5.9. In addition to assuming Hypothesis 5.2, suppose that p and x are 
periodic with the same period. Let D(p,x) = ( p -g- +x pS +X ) ^ e assocm t e d 
self-adjoint Dirac difference operator in l 2 (%) m ®i 2 {j€) m and suppose that D(p,x) 
has uniform spectral multiplicity 2m. Then D(p,x) * s reflectionless and hence for 
all feeZ, 

T D (X, k) = h 2m for a.e. A G a ess (D(p, X ))- (5.83) 

In particular, assume that p and x are periodic with the same period, that D(p,x) 
has uniform spectral multiplicity 2m, and that cr(D(p,x)) = [— E+ 2 , — E 1 / 2 ] U 
[E 1 / 2 , E^ 2 } for some < i?_ < E + . Then p and x are of the form, 

p{k) = diag(pi(fc), . . .,p m (k)), 
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Pj(k) = ^ ( E + 2 - £ i E - /2 ) > 1 < j < m, fc G Z, (5.84) 
X(fc) = diag(xiO), . • -,Xm(k)), 

Xj(k) = \ « /2 + e^l /2 ) , 1 < j < m, k e Z, (5.85) 

w/iere 

ej-e{l,-l}, l<i<m. (5.86) 
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